Suppressing decoherence and improving entanglement by quantum-jump-based 

feedback control in two-level systems 
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We study the quantum-jump-based feedback control on the entanglement shared between two 
qubits with one of them subject to decoherence, while the other qubit is under the control. This 
situation is very relevant to a quantum system consisting of nuclear and electron spins in solid states. 
The possibility to prolong the coherence time of the dissipative qubit is also explored. Numerical 
simulations show that the quantum-jump-based feedback control can improve the entanglement 
between the qubits and prolong the coherence time for the qubit subject directly to decoherence. 

PACS numbers: 73.40.Gk, 03.65.Ud, 42.50.Pq 



INTRODUCTION 



Superposition of states and entanglement make quan- 
tum information processing much different from its classi- 
cal counterpart. But a quantum state would unavoidably 
interact with its environment, resulting in a degradation 
of coherence and entanglement. For example, sponta- 
neous emission in atomic qubits [1] would spoil the co- 
herence of quantum states and limit the entanglement 
time. 

Recent experimental advances have enabled individual 
systems to be monitored and manipulated at quantum 
level @. This makes the quantum feedback control re- 
alizable. Among the feedback controls, The homodync- 
mediated feedback 0, Q and quantum-jump-based feed- 
back controls have been proposed to generate steady 
state entanglement in a cavity @. These two feed- 
back schemes are Markovian, namely, a feedback infor- 
mation proportional to the quantum-jump detection is 
synchronously used. Besides, these control scheme can 
also be used to suppress decoherence [7l-[Toj. 

Meanwhile, researchers are looking for proper systems 
for experimental implementation of quantum information 
processing. Among the various candidates, solid-states 
quantum devices based on superconductors (Tlj and lat- 
eral quantum dots [l2j are promising ones, however, the 
decoherence from intrinsic noise originating from two- 
level fluctuators is hard to engineered [l3j]. For this rea- 
son, the nuclear spins have attracted considerable at- 
tention [l4[ due to its long coherence times [Hj]. But 
their weak interactions to others make the preparation, 
control, and detection on them difficult. Thanks to its 
intrinsic interactions with electron spins, electron spin 
can be used as an ancilla to access single nuclear spin. 
This naturally leads us to rise the following question: 
can feedback strategy be used to suppress decoherence, 
prepare and protect entanglement between the nuclear 
and electron spins by controlling the electron spin? In 
this paper, we will study this problem by considering a 
nuclear spin (as a qubit) coupled to electron spin (as the 
other qubit) that is exposed to its environment. We show 



that a Markovian feedback based on quantum-jump can 
be used to suppress decoherence, produce entanglement 
and protect it. 

The paper is organized as follows: In Sec. II, we de- 
scribe our model and present the dynamics in absence 
of feedback. In Sec. Ill, we introduce the quantum-jump- 
based feedback control and give the dynamical equation 
under the feedback control. The effect of feedback control 
on decoherence and entanglement is discussed in Sec. IV 
and Sec.V, respectively. Sec. VI concludes our results. 



II. MODEL 

Our system consists of a pair of two-level systems, 
called qubit 1 and qubit 2, where only the qubit 2 inter- 
acts with its environment. We present a scheme employ- 
ing quantum-jump-based feedback control on the qubit 
2 to affect the decoherence of the qubit 1 and increase 
entanglement between the two qubits. The Hamiltonian 
of the system reads 



hg(a+a 2 + a x a£). (1) 



The first two terms represent the free Hamiltonian of 
the two qubits, the last term describes their interactions 
under the rotating-wave approximation. u\ and w 2 are 
the transition frequency of the two qubits, respectively. 
g is the coupling strength of the two qubits. a z is the 
Pauli matrix, i.e., a z — |e)(e| — \g)(g\, and er + = |e)(g|, 
°- = \9)(e\. 

The state of this quantum system can be described 
by the density operator p which is obtained by tracing 
out the environment. The dynamics of open quantum 
systems can be described by quantum master equations. 
The most general form of master equation for the density 
operator is [lj| [13] 



P : 



:[H,p]+C(j>), 



(2) 



where H is the system Hamiltonian and £ is a super- 
operator defined by C(p) = E k j k (L k pLl - \h\h k p - 
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\pL\L)~), in which different k characterizes different dis- 
sipative channels. 

In our system, the first qubit is assumed to be isolated 
from environment. The decoherence comes from the 
spontaneous emission of the qubit 2 (the second qubit). 
This situation is of relevance to a system consisting of 
nuclear and electron spins in aforementioned solid state 
devices. The dynamics of such a system takes, 

i 11 

P= -r[H,p\+7(a^pa^ - -a^u^p- -pa+a 2 ). (3) 

Here af — I\ ® af . The second part of Eq. © describes 
the dissipation of our system with 7 the decay rate. 

Though the first qubit is assumed to be isolated from 
environment, it still loss coherence due to the coupling 
to the second qubit. The decoherence process can be 
showed by the decay of off-diagonal elements of the re- 
duced density matrix for the first qubit. 

In order to investigate this decoherence, we calculate 
the evolution of system density operator p and then trace 
out the second qubit to get the reduced matrix 




The diagonal elements are the populations in the excited 
and ground states of the first qubit. And the off-diagonal 
elements represent the coherence of the qubit 1. 



III. QUANTUM-JUMP-BASED FEEDBACK 
CONTROL 

Quantum feedback controls play an increasingly im- 
portant role in quantum information processing. It is 
widely used to create and stabilize entanglement as well 
as combat with decoherence 0, HI, [l(| ■ In our model, 
the second qubit is used as an ancilla through which the 
feedback can affect the dynamics of the first qubit, i.e., 
by employing a feedback control on the second qubit, we 
control the first qubit. The goal is to suppress the deco- 
herence of the first qubit and enhance the entanglement 
between the two qubits by a feedback control on the sec- 
ond qubit @. 

Our feedback control strategy is based on quantum- 
jump detection. The master equation with feedback can 
be derived from the general measurement theory In 
our paper, Eq.© is equivalent to 

P (t + dt)= £ n a (2>(t)fit(T). (5) 

a=0,l 

with 

fti(cft) = ^fjdb(J2 (6) 
n Q ^l-C-H+^a+^)dt. (7) 



When the measurement result is a = 1, a detection oc- 
curs, which causes a finite evolution in the system via 
ili(dt). This is called a quantum jump. Then the un- 
normalized density matrix becomes p a= i = a 2 p(t)a 2 dt. 
The feedback control is added by giving p a= \ a fi- 
nite unitary evolution, then p a =i become p a =i = 
Fo2p(t)a 2 dt. In the limit that the feedback acts im- 
mediately after a detection and in a very shot time (much 
smaller than the time scale of the system's evolution) , the 
master equation is Markovian, 

P = -TW,p]+7(F(72pa 2 h F i - a 2 p - -po^cr^ ).(8) 

Here F = e iH s and H f = -\H' f t } . We see that the 
operator Hf contains a relatively large operator H'^ mul- 
tiplied by a very short time tf (Markovian assumption), 
but the product represents a certain amount of evolu- 
tion, so it is convenient to discuss Hf instead of Hf and 
tf. Here Hf is a 2 x 2 hermit operator which can be de- 
composed by Pauli matrixes Hf = A x a x + A y a y + A z a z 
(A x , A y , A z are real numbers). So we have, 

F = h <g> e lA ' 3 = h® (cos \A\ + i^^-A ■ a). (9) 

\A\ 

Here a — (a x , c yi a z ) and A — {A x , A yi A z ) representing 
the amplitude of a x , a y and a z control. 

In order to understand the physical meaning of feed- 
back operator F, we rewrite it as F = l 1 ®e~ % ^ n ' cr where 
n = (sin 6 cos 4>, sinO sin (j>, cos 9) and a = (a x ,a y ,a z ), 
this feedback operator is equivalent to a time-evolution 
with evolution operator F = I\ (g) e lHf . And it is 
clear that the operator F rotate the Bloch vector of 
the second qubit with the angle u> around the ft axis. 
The relationship between the two forms of F are A x = 
— j sin 9 cos 0, A y = — % sin 9 sin <f>, A z = — jCos9. So a 
a x control {A y = 0, A z = 0) means rotating the Bloch 
vector with a certain amount of angle around the x axis 
of Bloch sphere, so does the A y and A z control. Different 
A represents different feedback evolution i.e., rotate the 
Bloch vector with a particular angle around a particular 
direction in the Bloch sphere. For simplicity, we discuss 
the a x ,a y , a z control one by one in the following. 

This control mechanism has the advantage of be- 
ing simple to apply in practice, since it does not need 
real time state estimation as Bayesian feedback control 
does[l8(. The emission of the second qubit is measured 
by a photo detector, whose signal provides the informa- 
tion to design the control F . In this kind of monitoring, 
the absence of signal predominates the dynamics and the 
control is triggered only after a detection click, i.e. a 
quantum jump, occurs. 

IV. DECOHERENCE SUPPRESSION 

Before investigating the influence of the feedback con- 
trol, we first analyze the evolution of our system with- 
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out control. Assume that the two qubits are initially in 
the same pure superposition state, for example, — 
75(l e )i + 7§(l e )2 + 19)2)- The corresponding den- 

sity matrix is, 



1111 

/'" = ''/< r = 7 I \ 1 1 1 
1111 



(10) 



We assign the Planck constant H to be 1, uii — u) 2 = oj in 
Eq.([T]), and g/uj = 1,7/w = 0.5. After numerical calcu- 
lation, we get the evolution of the density matrix for the 
first qubit without control. Since p eg — p* ge , Pee + Pgg = 1 , 
we only discuss coherence \p eg \ and excited state popu- 
lation p ee for simplicity. The evolution of \p eg \ and p ee 
without control is depicted in FigQ] (a) and (b) (dashed 
lines). 

In FigQ] (a), a fast decay of \p eg \ (dashed line) can be 
found. This demonstrates that the first qubit lost coher- 
ence due to the second qubit 's spontaneous emission and 
their interaction. Meanwhile, the first qubit lost energy 
due to couplings to the second qubit (FigfTJ (b) (dashed 
line)). The results also show that the populations in ex- 
cited state and ground state decay away. This is because 
the first qubit exchange energy with the second qubit, 
see Eq.(UJ. 

Now we add feedback control F to our system, the mas- 
ter equation then becomes Eq.©. Our system is initially 
in the state po, other parameters remain unchanged. 

We first analysis the a x control by choosing feed- 
back amplitude A x = ~ w, A y = A z = 
0. Note that when A y — A z — 0, the feed- 
back amplitude A x influence the system's evolution 
a period of 7r which comes from the term 
It can be analytically proved that 



with 

Fa-pa+F^ in Eq.© 



e tA x a X(T - p2(X + e 

and e %AyC,v o~ p 2 o " e 
and A 



i(.Ai+7r)<Ta; 



~p 2 cr + e 



i{A x +ir)a x 

— e i(A v +ir)<r v a - p 2(T + e i(A y +ir)a y 

under any A x and A y . Here p 2 is the reduced density 
matrix of the second qubit. The absolute value for the 
first qubit 's off-diagonal density matrix element evolves 
as showed in Figf2](a). The figure indicates that, for an 
appropriate feedback amplitude , A x 1.3 and A x w 1.9, 
the absolute value of off-diagonal element can be ev- 
idently enhanced compared with the uncontrolled case 
(A x = 0). That means the decoherence is partially sup- 
pressed. The improvement of coherence caused by feed- 
back is shown explicitly in FigQ] (a). We plot \p eg \, rep- 
resenting the coherence of the first qubit, as a function 
of time with A x — 1.2, A y = A z = (a selected con- 
trolled case). In comparison with the uncontrolled case, 
a stronger oscillation amplitude and longer dechoherence 
time appears. Meanwhile, the p ee decays slowly com- 
pared to the uncontrolled case as shown in Figf2](b). 

Similarly, the a y control is also able to slow down the 
decay of \p eg \- We make A y = ~ n,A x = A z = 0. 
The numerical results of \p eg \ is shown in Fig [2] (b). 
Unlike the a x and a y control, the a z control (A = 
~ 7r, A x = A y = 0) has no effect on the evolution 
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FIG. 1: (a) Time evolution of \p eg \ with and without con- 
trol, t is in the unit of — . The different curves correspond to 
A x = 1.2, A y = A z = 0(solid line) and A x = A y = A z = = 
(dashed line), for g/uj — 1,7/w = 0.5. The feedback control 
strategy results in an improvement in decoherence time evi- 
dently, (b) Excited state population p ee evolution with and 
without control for the same parameters with (a), the decay 
of excited state population is slower in the controlled scheme. 



of the system as shown in FigJ5] (c). 



e iA ' <r 'o--p2(T + e 



This is because 
The physics be- 
hind this result is that after emitting a photon, the con- 
trolled qubit must stay in the ground state with the Bloch 
vector pointing the bottom of the Bloch sphere, so the 
rotation around z axis does not change the Bloch vector, 
i.e., the state of the qubit remains unchanged. 

The present results show that decoherence of the first 
qubit can be suppressed by controlling its partner. The 
decoherence source in our system is the spontaneous 
emission of the second qubit, once the detector detects 
a photon, i.e. a quantum jump of the second qubit hap- 
pens, the feedback beam instantaneously act on the sec- 
ond qubit and then the first qubit is impacted through 
the coupling of the two qubit. The feedback control 
scheme can reduce the destructive effects of coherence 
and slow down the dissipation of energy. 

The control effect is relevant to the coupling strength 
g. When g is small, the first qubit is hard to be impacted 
by the second qubit, so it's hard to prepare, measure and 
control the state of the first qubit. As the interaction goes 
stronger, the effect of feedback control becomes more ev- 
ident. 

For the case discussed in FiglTJ the first qubit is dis- 
sipative. We found that when the control parameters is 
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FIG. 2: The evolution of absolute value of the first qubit's off- 
diagonal element with different control parameters, for g/uj = 
1, y/w = 0.5 and t in the unit of —. (a) The a x control for 
A x = ~ 7T, A y = A z = 0. (b) The cr H control for Ay = ~ 
7r, Ax = A z — 0. (c) The a z control for A z — ~ tt, A x — 
A y — 0. When the feedback amplitude is chosen to be about 
f .3 and 1.9 for both a x and a y control, the oscillation of off- 
diagonal element is remarkably enhanced. The a z control 
doesn't work in our model. 



A y = with the two qubits initially being prepared in 
the same states, the decoherence dynamics turns to a 
phase damping type. The population in ground state 
and excited state do not change, while the off-diagonal 
elements evolves in the same way as in the uncontrolled 
case. We show this in a bloch sphere[19] in FigfS] Here 
the reduced density matrix of the first qubit can be writ- 
ten by pi = h{I + P ■ a). We can get the polarization 
vector components P x — Tr(a x pi), P y = Ti(a y pi) and 
P z = Tr(a z pi). 




-1 -1 



FIG. 3: Polarization vector evolution in a bloch sphere 
for feedback amplitude A x = \,A X = A y = (solid line) 
and A y — A x — A y = (dashed line). The parame- 
ters are g/to — 1,7/w = 0.5, the initial state is = 
45(|e>i + | 3 )i)®^(|e> 2 + |< 7 ) 2 ). 



V. ENTANGLEMENT CONTROL 

Quantum feedback control has been recently used to 
improve the creation of steady state entanglement in 
open quantum systems. A highly entangled states of two 
qubits in a cavity can be produced with an appropri- 
ate selection of the feedback Hamiltonian and detection 
strategy 0, Ho[ • We will show that the quantum-jump- 
based feedback scheme can produce and improve entan- 
glement in our model. We choose the concurrence [2l| 
as a measure of entanglement. For a mixed state repre- 
sented by the density matrix p the "spin- flipped" density 
operator reads 

p = (a y <g> a y )p*(a y ® a y ) (11) 

where the * denotes complex conjugate of p in the bases 
of {\gg}, \ge), \ eg), |ee)}, and a y is the usual Pauli matrix. 
The concurrence of the density matrix p is defined as 

C(p) = max(VAi- ^f\~ 2 - \/% - >/a!,0). (12) 

where Xi are eigenvalues of matrix pp and sorted in de- 
creasing order Ai > A2 > A3 > A4. The range of concur- 
rence is from to 1, and C = 1 represents the maximum 
entanglement. 

In absence of spontaneous emission, i.e. 7 = 0, the 
system evolves without dissipation. We find that for the 
system initially in a separable state except \ip) — |e)i|e)2 
or = \g)i\g}2 (the eigenstates of system Hamiltonian 
H), an entangled state can be generated due to the in- 
teraction between the two qubits. The amount of entan- 
glement depends on initial states of the system and the 
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FIG. 4: (a) Conccurence as a function of time and A y . The 
system is initially in the state = |(?) 1 1 e) 2 , for the param- 
eters g/ui — 1,7/w = 0.5. (b) A controlled evolution for 



A 1/ 



0.57T, A x = A z = vs. uncontrolled case. The entan- 



glement is improved by choosing an appropriate feedback, t 
is in the unit of ^ for (a) and (b). 



coupling strength g. But when the spontaneous emission 
effect is taken into account, the performance of entangle- 
ment preparation get worse considerably. 

Now we investigate if our feedback control strategy can 
improve the entanglement preparation with the effect of 
spontaneous emission of the second qubit. The master 
equation with control is Eq.®. The effect of feedback 
control lies in different choices for the feedback param- 
eters A x , A y , A z , the coupling strength g and different 
initial states. Here we present two typical results with 
two different states. 

Our first choice is the initial state \tp) = | (7) 1 1 e) 2 with 
dy control for A y = ~ tt, A x — 0, A z — 0. The con- 
currence evolution is plotted as a function of time and 
feedback amplitude A y in FigQ] (a), and Fig0] (b) de- 
notes the concurrence evolution with a selected feedback 
amplitude compared with the uncontrolled case. We see 
that entangled states can be generated with any feed- 
back parameters, but it decreases with time because of 
the dissipative effect. When an appropriate feedback 
amplitude A y w 0.9 is chosen, the concurrence ampli- 
tude is remarkably enhanced, and the entanglement lasts 
for a long time. For the system initially in the state 
\ip) = |e)i|e)2 with a y control, the dynamics of the con- 
currence is shown in FigJS] (a). Note that in this case if 
there is no spontaneous effect, this is a steady state of 
the system, the density matrix elements does not change 




FIG. 5: (a) Conccurence as a function of time and A y . The 
system is initially in the state \ip) = |e)i|e}2, for the param- 
eters g/uu — 1,7/cj = 0.5. (b) The controlled conccurence 
evolution for A y = 1.2, A x — A z — vs. uncontrolled case, t 
is in the unit of — for (a) and (b). 



with time. FigfS] (a) demonstrates that the dissipation 
and feedback can produce entanglement. We show this 
explicitly in Figj5] (b) by choosing feedback amplitude 
A y = 1.2. We can see that for a proper feedback am- 
plitude, after an entanglement death, a larger amount 
entanglement is regenerated. 

The above results shows the feedback control strategy 
can be used to prepare and protect entanglement in our 
model. The effect of entanglement control strongly de- 
pends on the initial state. For a certain initial state, we 
found that the a x control and a y control has the similar 
effect but the a 7 control docs not work. 



VI. CONCLUSION AND REMARKS 

In this paper, we studied the effect of quantum jump 
based feedback control on a system consisting of two 
qubits where only one of them subject to decoherence. 
By numerical simulation, we found that it is possible to 
suppress decoherence of the first qubits by a local control 
on the second qubits. We observed that the decoherence 
time of the first qubit is increased remarkably. The con- 
trol scheme can also used to protect the entanglement be- 
tween the two qubits. These features can be understood 
as that the feedback control changes the dissipative dy- 
namics of the system through the quantum-jump opera- 
tors. We would like to note that Hamiltonian Eq.(TT]) does 
not describe the hyperfine interaction. However, by the 
recent technology we can simulate Hamiltonian Eq. (TTJ) in 
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nuclear-electron spin systems, in this sense, the scheme tion Hamiltonian, our further simulations show that we 
presented here is available for nuclear-electron spin sys- can obtain results similar to that for Hamiltonian Eq. (JT|) . 
terns. On the other hand, by using the hyperfine interac- 
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